Spin entangled two-particle dark state in quantum transport through coupled 

quantum dots 
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We present a transport setup of coupled quantum dots that enables the creation of spatially 
separated spin-entangled two-electron dark states. We prove the existence of an entangled transport 
dark state by investigating the system Hamiltonian without coupling to the electronic reservoirs. 
In the transport regime the entangled dark state which corresponds to a singlet has a strongly 
enhanced Fano factor compared to the dark state which corresponds to a mixture of the triplet 
states. Furthermore we calculate the concurrence of the occupying electrons to show the degree of 
entanglement in the transport regime. 

PACS numbers: 03.65.Ud, 73.23.Hk, 73.63.Kv, 85.35.Ds 



I. INTRODUCTION 

The investigation of dark states (DS) has a 
long-standing tradition in quantum optics both 
experimentally 1 and theoretically 2 ' 3 . In recent years 
there have been numerous approaches to translate 
this quantum optical phenomenon into electronic 
transport 4-13 . Here, the term DS is used when the 
current-carrying particles, in general electrons, are 
trapped in a coherent superposition of states that is 
decoupled from the collector. The particle flow through 
the system is blocked, as no further electrons can enter 
the system due to the Coulomb blockade (CB). The 
first concepts used similar system setups as in quantum 
optics and included interactions with microwaves in 
order to create the DS 4,5 . 

A triple quantum dot (TQD) with a single excess elec- 
tron was the first system where an all- electronic DS was 
found by Michaelis et al. 7 ' 8 . Hence, the system is driven 
into the DS purely due to the coupling to the electronic 
reservoirs. Michaelis et al. showed the coherent trapping 
effect in the TQD and its destabilization due to charge 
fluctuations. This electronic DS was found to give rise 
to an enhanced Fano factor 8 above the Poissonian value 
F > 1. The influence of a magnetic field on the DS for- 
mation in the TQD was studied in Ref. [9,11] and Wey- 
mann et al. 14 have presented the effects of co-tunneling 
on the DS formation. The influence of phonon interac- 
tion on the dark state formation in the TQD was studied 
in Ref. [15], and Ref. [16] showed how the TQD dark 
state can be used as a nanomechanical resonator cooler. 

We have previously shown 17 that transport DS are not 
solely an issue of strong Coulomb blockade systems with 
only a single excess electron. In a TQD with a second ex- 
cess electron, a two-electron DS can be found for certain 
configurations. This two-electron DS can also be used as 
nanomechanical resonator cooler 18 . That electronic DSs 
also occur in interaction with other blockade phenomena 
is shown in Ref. [10] where a mixture of a spin blockade 
and a single-electron DS was shown to lead to a quasi 
two-electron DS. 



While the two-electron DS in the single TQD of 
Ref. [17] is a product state of two single-electron DSs, in 
this paper we introduce a system that enables the prepa- 
ration of a spin-entangled two-electron DS. For this aim, 
we consider two triple quantum dots with a single excess 
electron in each dot. A possible application of this setup 
is the creation of spacial separated entangled electrons 
on demand. 

The structure of this article is the following: After 
introducing the model in Sec. II, we investigate the exis- 
tence of dark states in the closed system without coupling 
to the electronic reservoirs in Sec. III. The transport 
properties, namely stationary current and Fano factor, 
are discussed in Sec. IV and in order to show the degree 
of entanglement in the transport regime we calculate the 
concurrence in Sec. V. 



II. MODEL 

Fig. 1 shows two possible configurations of the two 
TQDs. Both TQDs are in the strong Coulomb block- 
ade regime. Such that, up to one electron is allowed in 
each TQD. The TQDs are close together therefore we 
have a finite charging energy between the TQDs. Fur- 
thermore, we introduce an isotropic exchange interaction 
acting between the two TQDs. The complete closed sys- 
tem Hamiltonian Hd is given by 

Hd = -ffTQD.a + + U + J, (1) 

where the TQD Hamiltonians are given by 

3 2 

H T QB,A = ^2 2J Ei iA n lC r,A +Ta^2 *^2( d ]cr,A d 3v,A + h.c), 
i=l a j=l (J 

(2) 

with A <= {a, &}, difj A is the annihilation operator, and 
n-ia,A is the corresponding number operator of an elec- 
tron in quantum dot (QD) i of TQD A with spin a. We 
assume spin-independent energy levels and denote the en- 
ergy of the single-electron level of a quantum dot QDLA 
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(a) Triangular TQDs above each other (b) Serial TQDs in a parallel configuration 



Figure 1: Both TQDs are in the strong Coulomb blockade regime such that up to one electron is allowed in each TQD. The 
electrons in the TQDs interact with each other capacitively due the charging energies f/y and due to the exchange interaction 
switching the spins of the electrons. (The exchange interaction Ji is indicated in the sketches). Each TQD is connected to two 
sources and one drain. Here, two possible configurations are shown (a) where the two TQDs are triangular and lie above each 
other and (b) where the TQDs are serial and parallel to each other. In configuration (a) the special distance between QD1 and 
QD2 is smaller than in configuration (b), therefore the effects of decoherence which may destroy the DS in this configuration 
are smaller. 



as Ei.A- In the following we set all Ei A = 0. The levels 
in QD1,A and QD2,A are coupled coherently to QD3,A 
with a tunnel amplitude T A . The charging energy is the 
capacitive part of the Coulomb interaction and given by 

3 

U = '^2*^2Ui 3 h iai ahj r7 > i b , (3) 

where fZy is the additional charging energy needed to 
add an electron to QD i of TQD a when QD j of TQD 
b is occupied with one electron. In this setup an elec- 
tron in TQD a always interacts with an electron in TQD 
b, as both TQDs are in the strong Coulomb blockade. 
Therefore terms for having two electron in a single TQD 
are not included U, as they are assumed to far above the 
transport window and not relevant for the transport. 
The isotropic exchange energy is 

■I X ■ a 3,b)=Ys Ji i( a Z" a h +a ia a lb+°la a j,b) 

= J V ( a la a lb + 2 (°i, a a j,b + a i,a a t,b)) 

i,3 

= ^2 Jij [(flif, a ~ fki.,a){fkt,b - flii,b) 

*,j 

+ 2 (4 T ,o d ^.« d iW> d Jt.«» + 4i,a d ^,ad] f b d 3Ub )] , (4) 

where the er are the Pauli-matrices, i and j label the QDs 
of the TQD a or b and are the exchange constants. 
In the following we set Ja = J\ and = J2, i =/= j. In 



this article we treat this exchange interaction as a part 
of the Coulomb interaction between the electrons. 19 The 
system Hamiltonian in the localized basis can be found 
in Appendix A. We will later on refer to the two-electron 
Hamiltonian blocks H aa o'b for the different spin configu- 
ration and the exchange interaction blocks J defined in 
this Appendix. 

Each TQD is connected to three electron reservoirs 
that are described with the Hamiltonian 

Hres = ^ ^ e a k,AC aka A C a ka,A, (5) 
a, A k,a 

where a := {1,2,3} labels the reservoirs (1,2 = source, 
3 = drain) and <^ aka ^ is the creation operator of an elec- 
tron with spin a in mode k of reservoir a of TQD A. 
The TQD and the reservoirs are connected by the tunnel 
Hamiltonian 

Ht=J2Y1 v ak.Aci k ^ A d a<T , A + he (6) 

We assume spin-independent reservoir energies eik,A and 
tunneling amplitudes Vik,A- 

III. CLOSED SYSTEM 

The formation of a transport DS \&d) is only possible 
if the system Hamiltonian H$ fulfills certain conditions. 
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But the formation of the DS can be destroyed due to 
decoherence 4,20 or avoided for special coherent system- 
bath couplings 20 , even when these conditions are fulfilled. 
In the high bias regime 21 ' 22 all relevant system states 
lie well within the transport window. Here, a transport 
DS can be found when the system Hamiltonian block 
with most excess electrons has an eigenstate \&d) without 
finite occupation on the QD(s) which is (are) coupled to 
the collector(s). For the two-TQD setup this means, we 
search for an eigenstate in the two electron sector without 
occupation on QD3 of both TQDs 

(^MN = 0- VA,f. (?) 

In the transport regime such an eigenstate leads in gen- 
eral to a current blockade, where the stationary current 
of the system drops to zero when the DS becomes occu- 
pied. 

A. Closed system without exchange interaction 

We begin by looking at the system without exchange 
interaction and set J\ = J2 = 0, such that the charging 
energy is the only influence that exists between the two 
TQDs. Without exchange energy the two-electron sector 
of the system Hamiltonian consists of four blocks which 
are uncoupled to each other. Each block corresponds to 
one of the four possible spin configurations of the occu- 
pying electrons. In order to find a transport DS, we are 
searching for eigenstates without occupation on the third 
dots in one of these four blocks, hence eigenstates of the 
form 

\*D) =(aidL,a4a',b + a Aa, a dL',b 

+ MLa<<,6 + ML^L^IO), ( 8 ) 

with |ai| 2 + |a 2 | 2 + |a 3 | 2 + |a 4 | 2 = 1. But these blocks are 
9x9 matrices and it is in general not possible to calculate 
all eigenstate analytically, in order to prove that a dark 
state exists. However, the state \&d) must fulfill in the 
spin-erer' sector of the localized basis the condition 

(H <ra * b -\ D 1)\9 D ) = 0. (9) 

The explicit form of the H aa a'b can be found in Appendix 
A. Each of these four blocks has equal entries, because we 
assume spin degenerate single particle energies. Hence, 
a DS in one of the blocks is degenerated with the DSs at 
the same energy A^ in the other three blocks. 23 In the 
two-TQD setup without exchange interaction we find a 
DS 

law) = \(A a ,a 4,J(<4',6 - 4^)10), (10) 

in each block, when the charging energies U12 = Un = 
U22 are equal to the DS-eigenenergy Ad = Ui2- The 
DS of the transport system is then a mixture of the four 
degenerate states of the closed system. 



B. Closed system with exchange interaction 

With a finite isotropic exchange interaction the two- 
electron blocks with opposite spins couple to each other, 
while the blocks with equal spins remain uncoupled. Now 
a two-electron DS is found either when 

(H aai7b + J- \ D 1)\V D ) =0, or 

with 

\*d) =( a l dt l t ,a d H,6 + a *4.t,a4l,b + a 3 d h,a d H,b 
+ a i d h,a4i,b + b l d U,a d lt,b + b Al,a d h,b 
+ 63< iO 4t,6 + ^<,„4t,6)|0>. (12) 

and M 2 + M 2 + M 2 + M 2 + \h\ 2 + \b 2 \ 2 + N 2 + 
|6 4 | 2 = 1. DSs of the form of Eq. (10) still exist, when 
the two occupying electrons have equal spin, but now 
the charging energies have to fulfill the condition Un = 
U22 = J2 — Ji + U12 and the eigenenergy is shifted to 
X D = J2 + Ui 2 . 

In the opposite spin sector we find now two DS with 
different condition for the charging energies and with 
different eigenenergies. For \n = —3J2 + U12 and 
Un = U22 = 3 Ji — 3J2 + U12 the DS is a singlet state 

\*D,-) = -j=(\¥D,n)-\#DM)), (13) 

and for Xd = J2 + U12 and Un = U22 = J2 - Ji + U12 
the DS is a triplet state 

Wd,+) = \(\*D,n) + l#iut»- ( 14 ) 

These two dark states are entangled with respect to the 
spin of the electrons. In the transport regime only the 
singlet-DS can be prepared as a pure state, since it is not 
degenerate with any other DSs. The entangled triple-DS 
Eq. (14) is degenerate with the two DSs in the equal spin 
sectors, which also correspond to the other two states of 
the triplet. Not only the degeneracy of the eigenstate is 
broken, but also the conditions for charging and exchange 
energies is different for singlet-DS and triplet-DS. This 
enables to prepare the pure transport singlet-DS in the 
high bias regime, where all relevant system states are well 
within the transport window. If the degeneracy is bro- 
ken, but the condition for the charging energies remains 
equal for all four DS, it would be only for certain finite 
transport window configurations possible to prepare the 
singlet-DS as a pure state. 

IV. TRANSPORT PROPERTIES 

The considered transport setup is such that all elec- 
trons enter the TQDs from the source leads with the 
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(c) Steady state current for J=0ATc 



(d) Fano factor for J=0ATc 



Figure 2: (a) Steady state current (J) /Y as a function of the exchange energy J and SU normalized by Tc for Y/Tc = 1. 
The thin dark line (zero current line) in the density plot indicates the formation of the singlet-DS and the thick dark line 
the formation of the triplet-DS. Fano factor (b) (corresponding to (a)) is highly super-Poissonian around both DSs. But the 
Fano factor near the singlet-DS is around four times higher than around the triplet-DS. (c) Current as a function of SU/Tc, 
for different coupling strength Tc /Y at J = QATc- The curve for Tc = YY corresponds to a section through the density plot 
above. The Fano factors (d) show that the maximum value of the Fano factor near both DS is independent of the ratio Tc/Y. 
Parameter: V = 10Tb. 



rates TiA, i € {1, 2}, A € {a, b} depending on the QD and 
TQD in which the electrons tunnel and leave the TQDs 
by tunneling into the drain lead with the rate T^a- We 
assume that all considered energy levels of the two TQDs 
lie well within the transport window. We can therefore 
use a generalized master equation in Lindblad form 21 ' 22 
to described the transport through the two-TQD setup 



x \ 



;P D x D 



>x 



(15) 



The explicit form of the 12 coupling terms X £ ^\ Aa, 
with j G {1, 2, 3}, A € {a, b}, a e {t)4J; can be found in 
Appendix B. In order to calculate stationary current and 
the second-order zero-frequency Fano factor we rewrite 
Eq. (15) in Liouville space p{\) = (Wo + Je lx )p{x) an d 
introduce a counting field 24-27 . This counting field en- 
ables one to introduce the cummulant generating func- 
tion of the current distribution 



o) p(to)} 



(16) 



where Ttd{' • • } corresponds to the trace of the den- 
sity matrix. The n-th order zero frequency current 
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correlation 27 is then evaluated by 
d d n 



dt d{ixY 



Hx,t)\x=o,t- 



(17) 



The second order zero-frequency Fano factor is then de- 
fined as 



F = 



2)\ 



(Si 1 



(18) 



the 2-th order current correlation functions, normalized 
by the stationary current ((S' 1 )) = (J)). Since we 
are interested in the total current and noise through 
the system, we count the electrons tunneling from both 
TQDs. With this the jump operator becomes J p = 
a D3AapD\ Aty and Wo correspond to the other terms 
of the Lindblad equation Eq. (15). We could also count 
the electrons leaving each TQD separately. However, for 
the parameter setting, for which we calculate the steady 
state current and Fano factor in this article, the results 
would be simply half of the total current and Fano factor. 

In the following discussion of steady state current and 
Fano factor we set J\ = J, J 2 = 0, T a = Xf, = Tc, 
Tia = T 2A = T 3A = T, U u = U and U t] = V for 
i j. We then introduce SU = U — V as difference 
between intra-charging and inter-charging energy. Fig- 
ure 2 shows total steady state current and Fano factor 
of the two-TQD setup. Fig. 2(a) is a density plot of 
the current (I)/T as a function of exchange interaction 
J and charging energy difference SU normalized by Tc- 
The formation of both dark states singlet-DS as well as 
triplet-DS can be seen as dark lines running through the 
density plot. However, the width of the anti-resonance 
in the current around the triplet-DS is for this param- 
eter setting much broader than for the singlet-DS. At 
SU = J = 0, where the current valleys cross each other, 
the singlet-DS and the triplet-DSs of the closed system 
as well as of the transport system live in a degenerated 
subspace. The transport DS at SU = J = is therefore 
a mixture of the three triplet-DS and the singlet-DS. 28 
Fig. 2(b) shows the according Fano factor to the cur- 
rent density plot. Although the Fano factor reaches near 
both dark states highly super- Poissonian values, the max- 
imum value around the singlet-DS is strongly increased 
compared to triplet-DS. 

Fig. 2(c) shows the current as a function of the charg- 
ing energy difference SU/T, for different ratios of Tc/T. 
The current increases asymptotic with increasing Tc- 
Hence, the not shown current and Fano factor curves for 
Tc = 100r almost coincide with current and Fano fac- 
tor curves for Tc = 10T. The width of the current valley 
around the DSs decreases for both DS with increasing Tc, 
but again remains finite for Tc — > oo. Apart from that, 
the valley around triplet-DS decreases stronger than the 
valley around singlet-DS. 

In Fig. 2(d) the corresponding Fano factor to Fig. 2(c) 
is shown. The maximum value of the Fano factor near 
both DSs is independent of the ratio Tc/T and highly 
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Figure 3: (a) Concurrence as function of the exchange energy 
J and SU normalized by Tc for Tc = IF. (b) Concurrence as 
a function of SU normalized by Tc at J = OATc for different 
ratios of T c /r. Parameter: V = 10Tb. 



super-Poissonian. But the value of the Fano factor 
around singlet-DS is around four times higher than the 
value around triplet-DS. Similar features are found by 
Burkard et al. in 29 . The width of the Fano factor res- 
onance is widest for Tc = 1/2T and decreases for both 
smaller and higher values of Tc/T, in the shown plot. 
Therefore, the width of the Fano factor resonance is not 
simply decreasing for higher values of Tc/T as for the 
corresponding current. 



V. CONCURRENCE 

Entanglement is a very important aspect of quantum 
mechanics. It is responsible for the non-locality of quan- 
tum mechanics, which can be tested 30,31 via the viola- 
tion of the Bell's inequality 32 . Apart from that, coher- 
ences are the foundation of various concepts in quantum 
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mechanics such as quantum computation 33,3 , quantum 
teleportation 34-36 and quantum cryptography 34 ' 37 . 

A way to measure the entanglement of a mixed state 
is to calculate its concurrence 3 ® C . The concurrence 
of a two qubit system is define as C = max[0, yCVj" — 
y^j—o \f\j], with Xj being the eigenvalues of P2q{o~ v ® 
CTj / )P2q(c j , (8) o~ y ) in decreasing order. Here, p2Q is the 
density matrix of the two qubit system in the localized 
basis and the a y are the according Pauli matrices of the 
qubits. The spin degree of freedom is entangled in the 
two-TQD setup. In order to calculate the concurrence 
with respect to the spin qubits, it is necessary to trace 
out the QD states of the stationary state p s t a t of Eq. (15) 
Pspin = TrQ D [p st at]- The two particle sector of p spin cor- 
responds then to P2Q- 39 

Fig. 3(a) shows the concurrence as density plot for the 
same parameters as 2(a). At the pure entangled singlet- 
DS the concurrence rises to one indicating a maximal 
entangled states. Around the mixed triplet-DS the con- 
currence is zero which corresponds to a state without 
entanglement. Fig. 3(b) shows the concurrence as a func- 
tion of SU/Tc at J = 0.4T C for different ratios of T c /T. 
As for the according Fano factor, Fig. 2(d), the concur- 
rence resonance is widest for Tc = hT. 

VI. CONCLUSIONS 

We have shown that the preparation of a spin- 
entangled two-electron DS is possible. For this aim we 
have introduced a setup with two TQDs. The isotropic 
spin-exchange interaction between the occupying elec- 
trons lifts the degeneracy of the singlet-DS with the 
three triplet-DSs and enables the creation of a pure spin- 



entangled singlet-DS. The concurrence, which rises to 
unity at the singlet-DS, proves the existence of the entan- 
glement in the transport regime. Futhermore, the singlet- 
DS has a strongly enhanced Fano factor compared to the 
triplet-DS. This signature of the singlet-DS can be used 
to separate the entangled DS from the non-entangled-DS 
by measuring the Fano factor. 

As the electrons are still localized in the TQDs, this 
setup enables the creation of spacially separated spin- 
entangled electrons on demand. Without any futher 
modifications of the device the entangled electrons are 
simply stored in the two-TQD setup. Switching the 
chemical potential of the sources leads, such that they 
also become collectors, enables the usage the entangled 
electrons outside of the device. Here, the disadvantage 
is that each electron has two possibilities to tunnel out 
of the device, namely the two former source leads of 
the TQD it is occupying. Theoretical this can be eas- 
ily avoided by switching the tunnel rate of one of the 
former sources of each TQD to zero e.g. T2a = ^2b = 0. 
Experimentally it would probably be easier to consider 
a setup, which has only one source and one drain lead 
for both bias configurations. This does not change the 
essential features as the DS formation and the values of 
Fano factor and concurrence around the DSs. 
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Appendix A: Hamiltonian 



The system Hamiltonian of the two TQDs above each other has a block structure 
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where the zero in the first diagonal entry denotes the empty state, the next four entries of the form H a A denote the 
single-particle sectors of the two-TQD system, with A <E {a, b} labeling the TQD and a € {f, 1} labeling the spin of 
the electron. In the basis {|1a<t), ^Ait), |3a<t)} these part have the form 




(A2) 
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Here, Ta is the coupling term between QD 1 and 2 to QD 3 and 2Aa a detuning between the first and the second dot 
Ei.a = ^-Ai E\ a = —Aa and E$ t A = 0. The last four diagonal terms H aaa n, are two-particle sectors. In the basis 



the Hamiltonians become 

/ E/n-t-£A a -t-iAb U Jfc U U U U l a U \ 

/ C/i2+A a -Ai, T b T Q \ 



The J denotes isotropic exchange energy terms, with 
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The off-diagonal J terms in the Block structure Hamiltonian switch the spin of the electrons and the diagonal terms 
change the effective charging energy of the sector. 



Appendix B: Coupling terms 

We find 12 coupling terms Dx in the Lindblad equation Eq. (15) of the two-TQD setup: 



Dza<j — 


Vr 3aCT (|o)(3 Qa | 


ct' 


)(3aalba' \ ' 


f 2f )(T /)(3 a(T 2b -/ - 


- 3b (T /)(3 a(T 3b (T ' 


Dsba = 


Vr 36<r (|o)(3 6CT | 


-Ed 1 - 
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+ Ed^ 
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+Ed^ 
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- 3bcr')(2 a(T 3b (T ' 
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Vr 26CT (|o)(2 bCT | 


-Ed 1 - 


)(l aCT '2b (T - 


I ^aa' ) (^aer' ^bcr \ 





(Bl) 



In the following we set Tiao- = Ti.4, I^Ao- = r 2 ,4 and r 3j 4 CT = r 3j 4, A = a,b 1 a =T>4- We have assumed energy 
independent rates T^o- = 27r |Va fc(T: A| 2 <5(w - £ a k,A)- ^ 
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